Entangled cluster states of four interacting spins with frustrated quasiparticles 
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We constructed a four spin- 1/2 model on a triangular star lattice to study the quantum statistical 
behavior of the frustrated triple-spin clusters and double-spin clusters. The double-spin cluster 
behaves as boson in zero-energy state and the highest excited state. While in the ground state and 
other intermediate excited states, the double spin cluster shows a hybrid statistical factor of fermion, 
boson and exotic quasiparticles. We find a definition for computing the statistical factor between two 
triple-spin plaquette quasiparticles. The statistical matrix of exchanging two triple-spin clusters is 
far beyond fermion and boson. The eigenstates of this model are maximal entangled quantum states 
like Greenberger-Horne-Zeilinger state and Yeo-Chua's genuine four-qubit entangled state. The zero- 
energy state is a generalized singlet state of two double-spin clusters. While the highest excited state 
is generalized symmetric state of two double-spin clusters. The quantum entanglement between two 
double-spin clusters is relatively trivial in comparison with that of triple-spin cluster state. The 
von Neumann entropy of triple-spin cluster state is much higher. Majonana fermion has interesting 
physical behavior in strongly interacted fermion system. This spin model is not directly related 
to Majorana fermions. So we represent Majorana fermions by string operator of spins following 
the inverse Jordan- Wigner transformation. In this representation, the Majorana fermions do not 
behaves like independent free operators but entangled with one another. Exchanging two Majorana 
fermions not only changed the output wave function of the two fermions, but also rearranged all 
the rest components of the collective wave function. This property is showed by computing the 
statistical matrix of Majorana fermions for exchanging two particles. The output wave function 
of Majorana fermions is collectively entangled quantum state. The exotic statistical matrix maybe 
useful for constructing quantum gate. The current quantum information technology maybe can 
implement this four particle model. 

PACS numbers: 05.30.Pr, 03.67.Lx, 73.43.Cd 
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I. INTRODUCTION 

Entangled quantum states is the central source of 
quantum computation. Topological ordered states as 
a promising error-correcting quantum code for fault- 
tolerant quantum computation has aroused broad the- 
oretical research interest [H-[l3- There had been many 
proposals for experimentall y st udying quasiparticle in Ki- 
taev toric code model |l3|[l4[ as well as Kitaev honey- 
comb The laser excited Rydberg atoms in opti- 
cal lattice is an effective quantum simulator 11711 of Ki- 
taev toric code model and color code model (lq . The 
solid material Na2lrOs is also possible to map Kitaev 
honeycomb model into reality [19(. 

Majorana fermion appears as a promising candidate for 
topological quantum computation. It can be constructed 
by ordinary electron and hole in p-wave superconductor 
or quantum wire [20!] [21] . Direct observation of Majorana 
fermion is likely put into practice in quantum wire [22j . 
quantum do ts 1 231] , trapped ions |24j and vortex of atomic 
fermi gases |25j |. 

On the other route, quantum optics had been very 
successful to manipulate a small number of qubits. 
Greenberger-Horne-Zeilinger state has been generated in 
laboratory (26[, so did the W-state of four qubits [27| . 
Yeo-Chua's genuine four-qubit entangled state is also 
promising for experiment [28j . The current quantum 
information technology can implement these entangled 



eigenstates by cross-Kerr nonlinearity [29(. The excit- 
ing condensed matter systems above are somewhat too 
complicate for a quantum optical implementation. Thus 
we constructed a four-spin model with frustrated quasi- 
particles. It inherits the exotic physics of strongly corre- 
lated many body system, in the meantime it is convenient 
for quantum optical research. 

The strongly coupled four spins with frustrated quasi- 
particle is built on a triangular star lattice. The coupling 
rule between neighboring spins follows Kitaev honeycomb 
model. However we add an anti-ferromagnetic coupling 
between neighboring plaquette operators which sits right 
at the center of plaquette. This introduced the geomet- 
ric frustration of the three quasiparticles. A frustrated 
system has large degenerated Hilbert space and nontriv- 
ial topological order [3(j. Geometric frustration is the 
intrinsic reason for the exotic statistics of quasiparticles 
in this model. 

The article is organized as following: In section II, 
we solved the triangular star model by diagonalizing the 
Hamiltonian matrix. The eigenenergy values and spin 
configurations of eigenstates are derived. In section III, 
we computed the statistical factor of double-spin clusters 
and triple-spin clusters. In section IV, we analyzed the 
symmetry of zero energy state and the highest excited 
state and computed the von Neumann entropy of the 
triple-spin cluster state. In section V, we computed the 
permutation matrix of two Majorana fermion operators 
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on ground states. Exchanging two Majorana fermions 
goes far beyond fermionic statistics. The section VI is 
summary, we introduced the straight generalization of 
this triangular star model to a frustrated Kitaev honey- 
comb model. The appendix includes more details. 

II. THE ENTANGLED STATES OF TWO 
DOUBLE-SPIN CLUSTERS IN THE 
TRIANGULAR STAR MODEL 

In the triangular star model, the four particles are 
placed on the vertices of a triangular star lattice(Fig. 
[T]). The triangular star has three independent triangu- 
lar plaquette. The four particles coupled to each other 
following Kitaev honeycomb model. We added an anti- 
ferromagnetic coupling between the nearest neighboring 
plaquette on the Hamiltonian, 

+ JzG\G% + JpS\Si + J p S2S3 + JpSsS\. (1) 

The three plaquette operators commute with the Hamil- 




(a) (b) 

FIG. 1: (a) The triangular star model. Quasiparticle sits right 
at the center of each triangle, (b) The eigenenergy levels of the 
triangular star model with respect to the parameter setting 

Of (</p 2t/ a; , J Z 2t7a;, Jy St/a;). 

tonian. They are conserved quantum operators, 

Sx = a^al S 2 = ala y 2 a x 3 , S 3 = afojo*. (2) 

The three conserved plaquette operator divide the total 
Hilbert space into three sectors. Each plaquette opera- 
tor has eigenvalues +1 and —1 within its Hilbert space, 
S a \ip) = ±1|'0}. Every triangular plaquette operator de- 
fines an effective Ising spin within each sector. In the 
Kitaev honeycomb model, the ground state chooses a ho- 
mogeneous gauge pattern. All plaquette operators take 
the same eigenvalue. As there is antiferromagnetic cou- 
pling between two plaquette, the ground state is no longer 
the homogeneous gauge pattern. The ground state bear 
the frustrated gauge pattern. 

We first solve the model by diagonalizing the Hamil- 
tonian matrix. The spin operators take a sixteen dimen- 
sional representation, 

< = of ® I 2 ® I 3 ® I4, <?2 = !l ® & 2 ® r 3 ®l4, 

cr£ = Ii (g) I 2 ® a% <g> I 4 , (J% = Ii <g> I 2 ® I3 ® <t£, (3) 



where <rf are the conventional Pauli matrices and Ii is the 
2x2 identity matrix. The symbol <g> denotes direct prod- 
uct. The eigenvalues of the sixteen dimensional Hamil- 
tonian matrix lead to eight discrete energy levels, 

£± = 3J P ± 2yfjji + Jl + j], 

E± = -J p ±2J„, p = x,y,z. (4) 

Each energy level has two fold degeneracy. The energy 
levels are listed in Fig. [T](b). The eigenenergy and eigen- 
states are computed directly from the Hamiltonian ma- 
trix. The newly added coupling terms of two plaquette 
operators commute with Hamiltonian. It does not mod- 
ify the physics of Kitaev honeycomb model. The vor- 
tex excitation in the triangular plaquette represents the 
same type of quasiparticle excitation of Kitaev honey- 
comb model. The main different character of this tri- 
angle star model from Kitaev honeycomb model is the 
neighboring quasiparticles are now antiferromagnetically 
coupled to each other. So the ground state is the mini- 
mally frustrated quasiparticle state. 

The two levels with higher eigenenergy correspond 
to excited states for which the eigenvalues of the three 
plaquette operators takes the same value Si = S2 = 
S3 = ±1. The three pairs of energy level with lower 
energy correspond to the frustrated gauge pattern. 
Only one pair of quasiparticles is not frustrated for E^. 
If J x = J y = J z , the six energy levels of E^ would be 
degenerated. 

If the quasiparticle coupling interaction J p becomes 
zero, the triangular star model reduces to a finite Ki- 
taev honeycomb model on triangular star lattice. For a 
stronger quasiparticle coupling interaction than the spin 
coupling, 

Jp 2(/a;, J Z 2</a;, Jy 2 J x , (5) 

there exists four degenerated states with zero energy. 
Without losing important physics, we focus on this spe- 
cial parameter setting of Eq. ([5]) in the following. By rep- 
resenting the eigenvectors of the 16 dimensional Hamil- 
tonian matrix by the 16 four-spin basis, we derived the 
spin configurations corresponding to the four zero energy 
states, 
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(6) 



Here the symbols in the quantum wave function represent 
the double spin configurations, 

I ft) = Itt), 1 4) = 111), |o) = in), |o) = ut). (7) 

The eigenenergy of the states above is zero, i.e., H\o)^ = 
0, fj, = 1, 2, 3, 4. The zero energy state pair of \o}2 and 
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0)3 comes from the minimal frustrated gauge pattern, 
i.e, E£ = —J p + 2J X = 0. The two zero energy states 
of 1 0)1 and 1 0)4 correspond to the fully frustrated states 
E- = 3J P - 6J X = 0. 

The ground state has four fold degeneracy. We denote 
them as \g)y., fi = 1,2,3,4., 



\g)i = 

lff>2 = 
1.9)3 = 

\g)i = 
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(8) 



The eigenenergy of the four ground states are E g = 
-6J X , i.e., H\g)„ = -6J x \g) M V = 1,2,3,4. The 
ground states are the minimally frustrated quasiparti- 



cle states. Two of them comes from E y = — J p 



2J y — 



J rj 2J Z 



—6J X . The other two states is for E~ 
—6J X . The spin configuration of the two states \g)i and 
\g)3 is the superposition of three-up-one-down state and 
three-down-one-up state. The two states of \g)% and | (7)4 
is the superposition state of four spin basis with two-up 
spins and two-down spins. If we flip the spins of all of 
the four ground states, it generates a minus sign upon 
the original wave function. So the ground state breaks 
Z2 symmetry. 

The whole Hilbert space of eigenstates can be classified 
into two classes. The first class breaks the Z2 symmetry. 
This class only includes the ground state and the four 
degenerated states with respect to 2J X . The other class 
keeps Z2 symmetry. This class includes the zero energy 
states, the highest excited states and all the rest excited 
states. 

The energy level u = 12 J x is the highest energy level 
corresponding to the fully frustrated quasiparticle state. 
The spin configuration corresponding to this two-fold de- 
generated level is |e) 15 and \e)\e, 



|e)i 5 = |frfr) + l44)+2|00)+2|00) 

+ I 00) + I 00), 

|e)i B - \ U} + \ 44)+2|00)+2|00) 

+ 5| 00) + 5| 00), 



(9) 



The eigenstates of this model can be viewed as en- 
tangled quantum states of two double-spin clusters. For 
example, the nearest excited state above the zero energy 
state has four-fold degeneracy. The spin configuration of 
the four states are 



|e)ii = -|44) + |44), 

|e)i2 = -| ■frO) - I fiO) + I 40) 

|e)i3 = -| Of) + 1 04) - 1 o-fr) 

\e) u = -I 00) + I 00). 



40), 
04), 



(10) 



The corresponding eigenenergy with respect to {|e)n, 
|e) 12 , |e) 13, |e)i4} is E = 2J X . The spin configurations of 
|e)i2 and |e) 13 is the superposition of three- up-one-down 
and three-down-one-up. The state |e)n can be viewed 



as a dual state of the well-known Greenberger-Hornc- 
Zeilinger state for four qubit 26] , 



\GHZ) = \ff) + \ U). 



(11) 



|e)i4 is a singlet pair state of two double-spin clusters . 
The other two states, |e) 12 and |e)i3, can be viewed as 
generalized W-state of four qubits |27| . 



\W) = (| fO) + I fO) + I Of) + 1 Of)). 



(12) 



The two degenerated states of energy level E = —4J X 
are |e)g} and {|e)io, they both bear simmilar state struc- 
ture as W-state, 



|e) 9 = I fO) -I Of) -I 04) + |^0), 
|e)io = I fO) -I Of) -I 04) + 140). 



(13) 



These eigenstates are genuine entangled states of four 
spins. The zero energy states and highest excited states, 
{\o)i, 10)4} and {|e)i5, |e)i6}, can be classified into the 
same class as Yeo-Chua's genuine four-qubit entangled 
state (Hi, 



lx 00 ) - I 44) - I 4ft) - I 00) + I 00) 
+ |oo) + |oo) + |fr4) + lfrfr). 



(14) 



The eigenstates of this model can be implemented by 
the same operation as that for generating Yeo-Chua's 
genuine four-qubit entangled state using cross-Kerr non- 
linearity [29(. These states can be mapped into graph 
states following the strategy of Ref. [3l| , the observable 
operators in graph state theory has a physical meaning in 
this quantum model. Usually Yeo-Chua's genuine four- 
qubit entangled state or W-state is constructed without 
a Hamiltonian, but here we can derive the eigenenergy 
of these states. One can decompose an arbitrary entan- 
gled state as the superposition of these eigenstates. The 
weight of each eigenstate is marked by its eigenvalue. 
This offers us a new angle to see the internal structure of 
the quantum entanglement. 

The Wootters's concurrence provide a convenient 
way to quantify quantum entanglement of four-qubit 
states 1321 . 



Tabcd = M*%®Z v b ®Z v c ®{r v d \r)f 



(15) 



This concurrence operator is happened to be a conserved 
quantum operator of this triangular star model. The 
product of any two plaquette operators is a string of four 
identical spin operators, 



S1S2 = o\o\a\o%^ S2S3 

q c _ -V -V -V „y 
D3D1 — a 1 cr 2 a 3 a 4 . 



'2"4"3' 



(16) 

The plaquette operators keep an arbitrary ground state 
vector within ground state. For example, the opera- 
tion of plaquette operator Si on the vector ground state, 



[91 



ls)i, 13)2 , b)3 , 15)4 ] , gives a matrix, 
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The operation of concurrence operator reads, f = S3 Si, 
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III. THE STATISTICAL PROPERTY OF 
DOUBLE-SPIN CLUSTER STATES AND 
PLAQUETTE QUASIPARTICLES 

As all know, switching two fermions would generate 
a negative sign in front of the collective wave function. 
What would happen if it switches two double-spin clus- 
ters? Usually a particle composed of two fermions be- 
haves like a boson. However some of the eigenstates of 
the triangular star model behaves like a fermionic wave 
function of two double-spin clusters. We define a permu- 
tation operator, P[12;34], 

P[12; 34]*(1234) = (3412) = 77* (1234). (19) 

If 77 = +1, the cluster state is bosonic state. If r\ = 
— 1, the cluster state is fermionic state. If 77 bear more 
complex structure, we might call the clusters as exotic 
quasiparticles. 

The four degenerated ground states include one 
femionic cluster state and one bosonic cluster states, 

P[12;34]| 5 ) 2 = (-l)|.g} 2 , P[12;34]| 5 ) 4 = (+1)|<7> 4 .(20) 

The rest two ground states can map into each other by 
the permutation operator, 

P[12;3A]\g) 1 = \g} 3 , P[12;34]| 5 } 3 = l^. (21) 

So we define a four dimensional vector of ground cluster 
state, 

*[ fl (i 2M) ] = M<7>1, l3>2 ,| 5 >3 ,\9)i] T , (22) 
then the permutation operates on the vector states, 



i , [12;34]tt [fl(12 
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(23) 



If we select two states out of the four ground states as 
basis for quantum computation, the statistical factor of 
^[si>S3] = I Is) 1 ' 15)3 ] T ^ s a P auu matrix rj = a x . For 
the vector wave function of &[ g2 , gi ] = [ 1.9)2, 1.9)4 ] T , 
the statistical factor of &[g 2 , gi ] is another Pauli matrix 
rj = tj z . The statistical property of the nearest excited 
state above zero energy is similar to ground state, 

P[12;34]|e)n = |e) n , P[12; 34]|e)i 4 = -|e) 14 . 
P[12;34]|e)i 2 = |e) 13 , P[12; 34]|e) 13 = |e) 12 . (24) 



The four degenerate zero energy states and the highest 
excited states are all bosonic cluster states, 



P[12;34]|o)i = \o) u z = l,2,3,4; 
P[12;34]|e),- = |e)„ j = 15, 16. 



(25) 



The first excited state has two degenerated states. Per- 
mutating two clusters produce a statistical factor beyond 
boson and fermions, 



P[12;34]|e> 9 
P[12;34]|e)i 



e l7r [afaf]| e ) 10 , 
e™[a%a x 4 ]\e) 9 , 



(26) 



For the vector wave function of the first excited states, 
^[eg.eio] — [ l e )9, l e )io V i the statistical factor of double 
spin clusters is 
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1 



(27) 



The three plaquette operator define three quasiparti- 
cles. These quasiparticles are eigen-excitations of the 
quantum model. To find out if the quasiparticles are 
fermion or boson, we define the same permutation oper- 
ator as above, 

P[S 1 ;S 2 ] = P[l;4]*(1234) = ^ iiaa] *(1234). (28) 

The output of the permutation operator P[Si; S2] on the 
four dimensional vector ground states *[ g (i23 4 )] is a non- 
trivial statistical matrix, 
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'[Si;S 2 ] 



-erf af 

-1 

erf -erf 

0-100 



(29) 



Following the same procedure, one may find the statis- 
tical matrix of P[5 2 ;5 3 ] = P[l; 3]*(1234), P[Si;S 3 ], 
and so on. If quasiparticlc Si complete a closed trajec- 
tory and returns to its starting point, then the statistical 
matrices satisfies the relationship, 

P[S 3 ; S 2 }P[S 2 ; Si]P[Si;S 3 ]P[Si;S 2 ] = I, (30) 

where / is a four dimensional identity matrix. 

In the first excited state, ^[e 9 ,e w ] — [ l e )9, l e )io ] T > the 
triple-spin cluster is a hybrid of fermion and boson. The 
fermionic part and bosonic part of the first state can be 
separated as 

A 9 = |K)) + |K>>, B 9 = -| Oft) - | OJ|), 
A w = -| Ofr) - | OJ1-), S10 = |1K5> + UO>. (31) 

The operation of permutation operator P[Si \ S 2 ] on the 
first excited states, |e)g = A 9 + P 9 , |e)io = A w + B w , 
gives us, 



P[5i;S a ]|e> 9 
P[5i;5 2 ]|e) 10 



-A10 + Bg, 
-Ag + BiQ, 



(32) 



■5 




FIG. 2: The permutating steps for plaquette quasiparticle Si 
to complete a circle around the center. 

We define the statistical factor of the vector wave func- 
tion of *[ e9)ei0 ] = [ |e) 9 , |e)i ] T as 

P[SilS 2 }^ [e9 , eio] =^ Si . s ^ [egteio] . (33) 
The statistical factor can be derived by matrix algebra, 
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(34) 



here J is 2 by 2 identical matrix. So the plaquette quasi- 
particle has complicate statistical behavior. This exotic 
property may be used to construct non-trivial quantum 
logical gate. 



Both |o)i and |e)i6 are bosonic state of double spin 
clusters. However the statistical factor of triple-spin clus- 
ter state is beyond fermion or boson statistics. Any 
two of the three plaquette always share two spins, so 
one can not cut one triple-spin cluster out of the four 
spins to get an independent three-spin wave function. 
We use the reduced density operator to represent the 
triple-spin cluster state. The density operator for the 
zero energy state |o)i of quasiparticle Si — <y 7 \02 a \ 
reads p Si (oi) 
lation yields 



Tr [ 4 ] [[I°)ii(°l][i234]]- A detail calcu- 
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We also derived the density operator of Si in the highest 
excited state |e)i6, 



IV. THE TRANSFORMATION BETWEEN 
SYMMETRIC CLUSTER STATE AND 
ANTISYMMETRIC CLUSTER STATE 

The conventional spin 1/2 particles may form a anti- 
symmetric state (| ti) — I it)) an d a symmetric state 
(I tt) + I ii))- Usually the singlet state leads to lower 
energy. The triplet leads to higher energy. A singlet can 
be transformed into a triplet state by flipping the second 
spin or the first spin, 



It) 



=> 



tt) + I ii) ]• (35) 



The singlet state and triplet state above is defined for 
two spins. While the spin configuration of the triangle 
star model includes four spins. A similar transformation 
rule as that for a conventional singlet state and a triplet 
state also exist between the zero energy state and the 
highest energy states of this triangular star model. Here 
is we need to flip a pair of spins. Based on the spin con- 
figuration of zero energy states and the highest excited 
states, one may extract two unit configuration out of the 
complete states, 



A = |M) + 2|00) + 2|00) + |M>, 
B = 5| 00) + 5| OO). 



(36) 



Both the zero energy state and the highest excited state 
are linear combination of the two units, 

\o)x=A-B, \e) 16 = [aM]A-e^[a^]B. (37) 

The same algebra relation also exist between the state 
\o)i and |e)i5 . The two unit spin configurations of [0)4 
and |e)i5 are different from that of |o)i and |e)ie- 
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(39) 



The two density operators of quasiparticle Si have the 
same eigenvalues, A = [1/2,1/2,0,0,0,0,0,0]. Only the 
two non-zero eigenvalues, Ai = 1/2 and A2 = 1, con- 
tribute to von Neumann entropy, 



S(p) = -AiLo0[Ai] - A 2 Lo ff [A 2 ] = 1. 



(40) 



von Neumann entropy quantifies the quantum entangle- 
ment of the triple-spin quasiparticle. The von Neumann 
entropy of quasiparticle S2 = tr|cr|(j| and S3 = afafa^ 
are the same as Si. The density operator of S2 for the 
zero energy state \o)i reads 
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Correspondingly the density operator for the highest ex- 
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cited states |e)i6 is 
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(42) 



The reduced density operator of S3 is originally based 
on the spatial order [124]. If we put the last spin No. [4] 
in front of [12], the modified spatial order of the three 
spins [412] generates the same density operator as that 
of 5*2. The eigenvalues of p s [124] (01) is the same as that 
Ps 3 [4i2] Thus the density operator of S3 is equivalent 
to 5*2, i.e., p s? (oi) = p S2 (oi), p S3 (e 16 ) = p S2 {e 16 ). 

Both the eigenvalues of S2 and S3 is the same as the 
quasiparticle Si . The von- Neumann entropy of the den- 
sity operator of the three quasiparticles are the same as 
well. Since the three quasiparticles intersect with one an- 
other, it is not a surprise if they have the same degree of 
quantum entanglement. 

V. THE STATISTICAL FACTOR OF 
MAJORANA FERMION OPERATORS ACTING 
ON GROUND STATE 

Kitaev honeycomb model was solved by two different 
approaches, the first one is representing a spin operator 
by two Majorana fermions [l[. The second approach is 
using Jordan- Wigner transformation Jordan- Wigner 
transformation maps spin to fermion or vice vers a[33j . 
Here we combine the two different approaches [7[ to 
represent Majorana fermion by spins. These Majorana 
fermion operators are composite operators, it generates 
complicate entangled states. 

The explicit formulation of inverse Jordan- Wigner 
transformation depends on the spatial order of the four 
particles. We squeeze the triangular star lattice to a one 
dimensional chain by keeping the topology of interact- 
ing bonds invariant. The indices ('1'/ 2',' 3'/ 4') in the 
Hamiltonian Eq. (jTJ) are the name of the four particles in- 
stead of its spatial ordering in the one dimensional chain. 
We denote the four spatial positions along the one dimen- 
sional chain as (ii < 12 < 13 < 14) (Fig. [3]). If we place 
particle '1' on 14, the other three particles would sit on 
the rest three sites (ii, 12, 13) before 14. For a special case 
of the spatial ordering of the four particles [1423], i.e., 
(ii = 1, '12 = 4, i3 = 2, 14 = 3), the inverse Jordan- 
Wigner transformation defines the spin representation of 
eight Majorana fermions, 

ipi = of, ip 4 = a\a\, tp 2 = (j\a\a\, 
h = -of, b 4 = -<jl<T z 1 , 6 2 = -ofo-iVf: 
■4)3 = ofofofo-f, b 3 = -of ojo-fo-f . (43) 




FIG. 3: The one dimensional scheme of the triangular star 
model. The initial spatial ordering of the four particle is 
[1423]. Then we braid the two particle at the second site 
and the third site, it leads to [1243]. The coupling bonds 
connecting the four particles remain the same as the original 
triangular star model. 

The string operator of spins for Majorana fermion exactly 
fulfills the commutator of the original Majorana fermions, 
b l = h, tpl = ipi, {bi,bj} = Sij, {ipi,ipj} = 5ij, {h,ipj} = 
0. 

First we need to label the spatial position of the Ma- 
jorana fermions. The inverse Jordan- Wigner transforma- 
tion provides a natural way for labeling each Majorana 
fermion. The first Majorana fermion is represented by 
one spin operator, the second is expressed by the prod- 
uct of two spin operators, and so forth. If we exchange 
the position of two Majorana fermions, its representa- 
tion would change spontaneously, so does its action on 
quantum state. 

Then we introduce a vector of eight identical |<?} 4 
states, 

\g)4 = [\9}i,\9)i,\g)^\9)i, [5)4! b)4 ; ls) 4 Js) 4 ] T ( 44 ) 

We define a spatially ordered matrix operator of the eight 
Majorana fermions 

FlhhisU] = diag[i)i,ip2,ipz,ip4,MMMM]- (45) 

The subscribe number T' of the fermion ipi is the name 
of particle, a t different position is represented by a 
string operator with different length. The subscribes of 
i 7 '[i 1 i 2 i 3 i 4 i defines the spatial order of the four particle. For 
example, ^4231] means particle No. '1' sits in the fourth 
site i4. Then ipi is represented by the product of four 
spin operators. This diagonal matrix operator always has 
the form diag[ipi, ip2, ip3, V^; b\, 62, &3, 64]- The position 
of a Majorana fermions is represented by the length of 
spin representation. When we exchange two Majorana 
fermions, the string representation of the two Majorana 
fermions would exchange. 

We operate the Majorana fermions on the ground state 
|.g) 4 to study their statistical matrix. The output vector 
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of matrix operator i'taiaijiJ on ^4 can be denoted as 



1% 



l%)2, |%>3 ! l%)4, 
l%>5, l%>6, l%>7, l%)8,] T - 



(46) 



The action of the matrix operator is equivalent to a sta- 
tistical matrix, 



(47) 



The output matrix of -f r [i 1 1 2 i 3 i 4 ] describes the statis- 
tical behavior of the Majorana fermions. If it is 
simply F = diag[l, 1, 1, 1, 1, 1, 1, 1] after we exchange 
two Majorana fermions, that means the two Majo- 
rana fermions behaves like boson. If it is F = 
diag[— 1, —1, —1, —1, —1, —1, —1, —1], that means the two 
Majorana fermions behaves like ordinary fermions. If 
-f]iii 2 i3i4] i s complicate matrix beyond the two above, the 
two Majorana fermions have exotic statistics. 

We set the initial position of the eight Majorna 
fermions as -F[i423] • Eq. (|43|) summarizes the spin repre- 
sentation of Majorana fermions for the spatial ordering 
of [1423]. We operate these string operators on \g) A and 
express the output spin configurations by linear combi- 
nation of other eigenstates, 



■y(l°>3 + |e>io), = ^(l«) 2 | + |e) 9 ), 

-y(|o)3 - |e)io), |%>4 = -y(|o)2 - |e) 9 ), 

|(|e>i2 ~ k%)e= 2(l2)i + |e)i 2 ), 

^(lff>3 - |e)i 3 ), l%)s = ^(|e)i 3 + \g) 3 ). 

(48) 



l%)3 

l%)s 



The eight dimensional output vector is the initial 

output for the spatial ordering [1423] . The initial matrix 
for F [1423 ] is an identity F [1423] = diag[l, 1,1,1,1,1,1,1] 
since the exchanging operation is still not performed yet. 

We take two different paths of permutation operation 
for the four particles. For the first path, we first exchange 
the second particle with the third particle, then exchange 
the third particle with the fourth particle. Here we in- 
troduce a permutation operator -Bi 2 i 3 for exchanging the 
particle at site i 2 with that at i 3 , i?i 2 i 3 [1423] = [1243], 



F 



[1243] 



/l 

















°\ 














—i 











-1 



































i 











1 











% 



































-1 





\0 





—i 











0/ 



(49) 



-PI1243] suggest that exchanging two Majorana femions 
does not only add a minus sign in front of the initial wave 



function. The we continue to perform the second permu- 
tation upon the first permutation, i?i 3 i 4 [1243] = [1234]. 
The final permutation matrix for the final ordering [1234] 
is 



F 



[1234] 



(I 









i 



\0 



-1 





—i 





i 



0\ 



1 





0/ 



(50) 



For the second path, we first exchange [1314] and then 
exchange [i 2 i 3 ] (Appendix E), 



£ i3i4 [1423] = [1432], £ i2i3 [1432] = [1342]. 
The final permutation matrix for [1342] reads 



(51) 



F 



[1 3 42] 



/I 




























—i 


























1 




























—i 




























1 




























—i 




























1 







Vo 




















—i 


) 



(52) 



Usually if we exchange Majorana fermions twice, the 
statistical matrix would be a positive identity matrix. 
However here the permutation matrix ^[12341 and -Fji 34 2] 
above are obviously non-trivial matrix comparing with 
identity matrix. In fact, these permutation matrix form 
a group, when the spatial order of the particles returns 
to its initial ordering after a series of permutations, the 
product of those permutation matrices along the path 
should be equal to identity. 

For a more general theoretical framework for permu- 
tating particles by Jordan- Wigner transformation, when 
a fermion moves one step forward, it would be multi- 
plied by a z . If a fermion moves one step backward, a a z 
would be taken away from its string operator. To explic- 
itly formulate this operation, we build ip^ fermion and 
fermion into a two dimensional vector, 



(53) 



There are two independent ways to switch the nearest 
neighboring two vectors of fermions operators. The first 
way is to map i/^-i to an d to 6s, the permu- 
tation matrix is 



V-i 

a. 







(54) 



The other way is transforming if^i _ x to &i M and mapping 
6i f4 _ 1 to tpi^, the permutation matrix is off-diagonal ma- 
trix 







V-i 




(55) 
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Both the two matrices satisfy the algebra relation, 

M^iii' ( 56 ) 

Every lattice site shoulders a pair of dual Majorana 
fermions. This is to compensate for the degree of free- 
dom of spin operators. The example above suggest that 
only when we put all of the eight Majorana fermions into 
action, the Hilbert space generated by the fermion oper- 
ators form a closed space. This suggest the eight Majo- 
rana fermion are highly entangled set of operators. This 
is another way to understand quantum entanglement. 

VI. SUMMARY 

The quantum statistics of exotic quasiparticles is more 
complicate than fermion and boson. It is a challenge to 
observe exotic quasiparticles in quantum optical experi- 
ment nowadays. We proposed a four particle model to 
study the explicit spin configuration of exotic quasiparti- 
cles which maybe is useful for quantum optical research. 

This triangular star model is inspired by Kitaev honey- 
comb model. We introduced the spin coupling rules in Ki- 
taev honeycomb model but add frustrated coupling terms 
between plaquette quasiparticles. Direct diagonalization 
of the Hamiltonian matrix produces the complete list of 
eigenenergy levels and eigenstates with explicit spin con- 
figuration. The eigenstates of this model are genuine en- 
tangled four qubit states. For instance, the first excited 
state can be viewed as a generalization of Greenberger- 
Horne-Zeilinger state and Yeo-Chua's genuine entangled 
states. 

The quantum entanglement and quantum statistics for 
a clusters made of several-spins depends on how we pack 
the spins. We first take double-spin as a cluster. Usually 
a composite particle made of two spin 1/2 behaves like 
boson. This is indeed the case in zero energy level and 
the highest energy level where exchanging two double- 
spin clusters keeps the collective wave function invariant. 
However when we perform the same exchanging opera- 
tion on ground state and intermediate excited state, there 
appears non-trivial off-diagonal elements in the statisti- 
cal matrix. A double-spin 1/2 cluster may be a boson, a 
fermion, and other exotic particle. The zero energy state 
may be viewed as a singlet state of two double-spins. The 
highest excited state is a symmetric state of two double- 
spins. 

The triple-spin cluster happened to be eigen- 
quasiparticles of this four spins model. We computed the 
statistical matrix for exchanging two triple-spin clusters 
on the first excited state, it does behave like a fermion as 
what we usually believe. The statistical matrix is a com- 
plicate matrix. So these triple-spin quasiparticles are nei- 
ther boson nor fermions, they are exotic quasiparticles. 
Even in the zero energy state and the highest excited 
state where double-spin are normal bosons, the triple- 
spin cluster is still highly nontrivial, the von Neumann 
entropy of the reduced density operator of the triangular 




FIG. 4: . The extension of frustrated plaquette excitation 
model to Kitaev honeycomb model. 

plaquette excitations is much higher than that of double- 
spin clusters. So the quantum entanglement in a triple- 
spin cluster is much stronger than that in a double spin 
cluster. 

Permutation operation is the elementary operation for 
constructing quantum logical gate. Permutating Majo- 
rana fermions is believed to be a promising approach to 
quantum computation. However a free Majorana fermion 
is absent in nature. So we construct a composite Majo- 
rana fermion by a string of spins in the inverse Jordan- 
Wigner transformation. It fulfils the commutator rela- 
tionship of Majorana fermions. The string operator of 
spins for Majorana fermion operator is different for dif- 
ferent spatial ordering. So exchanging two Majorana 
fermions is encoded in the spatial order of spin repre- 
sentation. When we exchange two Majorana fermions, 
the effective statistical matrix is not simply adding a mi- 
nus sign in front of the original wave funtion. It indi- 
cates the output states of a composite Majorana fermion 
is strongly entangled with other composite Majorana 
fermions. This statistical matrix is actually equivalent to 
permutation matrix. When the spatial order restored its 
original order after a circle of permutations, the product 
of these permutation matrices would be equal to identity. 
This operation constructs a permutation group. 

The frustrated coupling term is the source of those 
exotic statistical behaviors. It can be straightly ex- 
tended to Kitaev honeycomb model, we need to add anti- 
ferromagnetic coupling between the nearest neighboring 
plaquette upon the Kitaev honeycomb Hamiltonian Hk, 
H = Hk + Hi , 

H k = ^2(J x a^a^ +e3: + J y o v v a v r+t3y + J z a*a* +e J, 

r 

Hi = J2 Si, = trfofoMofolf, (57) 

<ipj P > 

where i p is the center of every plaquette(Fig. HJ. All 
plaquette operators commute with Hamiltonian as well 
as themselves, [S ip ,H] = 0, [S"i p , Sj p ] = 0. The newly 
added quasiparticle terms commute with Kitaev honey- 
comb Hamiltonian, [Hk,Hi] — 0. The plaquette coupling 
interaction is antifcrromagnetic, i.e., J p > 0. This model 
can be solved by inverse Jordan- Wigner transformation. 
A string operator of many spin operator is convenient 
for manipulating plaquette state. Further research will 
be carried out in the future. 
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VIII. APPENDIX 

A. The inverse Jordan- Wigner transformation for 
different spatial ordering 



iIj-]\q)a = 


— (|0>3 ■+ 


1 / J-uy 7 


02 1 Q)a 


= — (\g) 1 + 


|e)i 2 ) 


03I3U = 


— (I.9>3 - 


- |e)is), 


^4|5>4 




e)g), 


frilsU = 


>12- 


l.9>i), 


b2\g)i 


= ^(l°>2 + 


e> 9 ), 


b 3 \g)i = 


4d°)3- 


|e)io), 


bi\g)i 


= i\(\g)i + 


|e)i 2 ). 



(62) 



The inverse Jordan- Wigner transformation for the or- 
dering [1243] is 



The inverse Jordan- Wigner transformation for the or- 
dering [1324] is 



01 = erf, 

V 
T 2< 



b 3 = -ala{olal 
4 = — crfof af 

■04 = <yt<J z a 2 , b 2 = 

V-3 = crfofo-fof , 61 = -of. 



~ x ~z 

a 2 °1, 



(58) 



The output state of the eight fermions for [1243] on 
ground state | (7)4, 



0l|ff)4 
02|3)4 

04(3)4 

03|ff)4 



— % 1 

y (|o) 3 + |e)io), &ib>4 = ^ (I e >i2 - |3>i), 

y(|3>i + |e>i 2 ), b 2 \g)i = ^(|o) 2 + |e) 9 ), 

-|(l.9>3 + |e)i 3 ), b 4 \g) 4 = ^(\o) 2 -\e) 9 ), 

1 i 

2(|o)3 - |e)io), 6s|ff>4 = 2 d e )!3 _ lff)s)- 



(59) 



Comparing the output states above with the output state 
for the initial spatial ordering [1423], we derived the per- 
mutation matrix for the spatial ordering [1243], 



01 = crf, 64 = -a\alulal, 

03 = ofcf, & 2 = -erfef of , 

„/, „y „z „z 1 „y z 

02 = u 2 u x g 3 , b 3 = -a 3 a 1 , 

04 - olalalal, h = -of. (63) 

The output state of the eight fermions for [1324] on 
ground state \g)±, 



—i 1 

0i|3)4 = -y(|o)3 + |e)io), 03|3)4 = ^(\o) 3 - \e) 10 ), 

i 1 

02|#)4 = 2(10)2 + |e> 9 ), 4 |3)4 = ^(\o) 2 - |e) 9 ), 

1 i 
h i\g)i = 2(l e )i2 — bz\g)i = 2-(l e )i3 - b)s), 

02|5)4 = y (|.9)i + |e)i 2 ), 64|3)4 = y(|e>i3 + |ff>3)- 



(64) 



The inverse Jordan- Wigner transformation for the or- 
dering [1342] is 



'[1243] 



(I 

















°\ 














—i 











-1 



































i 











1 











i 



































-1 





\o 





—i 











0) 



(60) 



The inverse Jordan- Wigner transformation for the or- 
dering [1234] is 

01 = erf, 64 = -crfoferfof, 

02 = <y\o 7 i , 63 = -of of af , 

03 = Crf Of erf, 6 2 = -ofof , 

04 = of cf crfcrf , 61 = -erf. (61) 

The output state of the eight fermions for [1234] on 



01 = crf, b 2 = -a\a\Gla z 2 , 

03 = of of, b 4 = -of of crf, 

04 = crf cf erf, 63 = -crf crf, 

02 = a x 2 a\a z 3 a z 2 , b x = -crf. 



(65) 



The output state of the eight fermions for [1342]on 
ground state \g)i, 



-i -1 

0i|3)4 = ~j{\o) 3 + |e)io), 03|ff)4 = — (|o) 3 - |e)io), 

i 1 

04|3)4 = 2(10)2 - l e )g), 02|.9) 4 = 2(10)2 + |e) 9 ), 

1 —i 

bi\g)i = ^(\e)i2 - \g)i), b 3 \g) 4 = y(|e)i 3 - \g) 3 ), 

62|3)4 - y(l3)i + |e)i 2 ), & 4 |3)4 = y (|e)i 3 + |3)s)- 

(66) 
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B. The 16 basis of spin configuration 

The 16 basis for expanding the eigenvector of the 16 
dimensional Hamiltonian matrix, 



11) = 


tttt}, 


|2} = 


tut}, 


13} = 


tut}, 


14} = 


1 ttU}, 


15} = 


titt), 


|6} = 


Utt), 


17} = 


itut>, 


|8} = 


uu>, 


|9} = 


tut), 


|10) = 


um>< 


|H> = 


i tut), 


|12} = 


i tut), 


|13) = 


lUtt), 


|14} = 


i utt), 


|15} = 


1 tut}, 


|16) = 


1 -i-l^-l^) ■ 
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